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Abstract 

Neutrino-induced upward-going muons in MACRO have been anal¬ 
ysed in terms of relativity principles violating effects, keeping standard 
mass-induced atmospheric neutrino oscillations as the dominant source 
of —> Vt transitions. The data disfavor these exotic possibilities even 
at a sub-dominant level, and stringent 90% C.L. limits are placed on 
the Lorentz invariance violation parameter |Ati| < 6 x 10“^'^ at sin20„ 
= 0 and |Ari| < 2.5 5 x 10“^® at sin20i, = ±1. These limits can 

also be re-interpreted as upper bounds on the parameters describing 
violation of the Equivalence Principle. 


Neutrino mass-induced oscillations are the best explanation of the 
atmospheric neutrino problem [HElEllll. Two flavor oscil¬ 

lations are strongly favored over a wide range of alternative solutions 
such as Vsteriie oscillations ElEl, — > Ve oscillations 011 or 

other exotic possibilities 130 . 

In this letter, we assume standard mass-induced neutrino oscil¬ 
lations as the leading mechanism for flavor transitions and we treat 
Lorentz invariance flavor transitions as a sub-dominant effect 0. In 
the literature, neutrino oscillations induced by Violation of (CPT- 
conserving) Lorentz Invariance (VLI) and Violation of the Equivalence 
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Principle (VEP) are described within the same formalism. In the fol¬ 
lowing we will mention only VLI for simplicity. 

In this scenario, neutrinos can be described in terms of three dis¬ 
tinct bases: flavor eigenstates, mass eigenstates and velocity eigen¬ 
states, the latter being characterized by different Maximum Attainable 
Velocities (MAVs) in the limit of infinite momentum. 

Both mass-induced oscillations and VLI transitions are treated in 
the two-family approximation and we assume that mass and velocity 
mixings occur inside the same families (e.g. \v 2 ) and It's)). 


The usual interpretation of the atmospheric neutrino oscillations is 
Vt induced by the mixing of the two mass eigenstates liz™) and 
\v^), and two weak eigenstates \Vfj) and \vt)^ he. 

= \vir) cos 9^3 + WT) sin^ 2 ! 

Wr) = -K) sin^al + WT) cos0^3 

where 0^ (= 9m) is the flavor-mass mixing angle. The survival prob¬ 
ability of muon neutrinos at a distance L from production is 


P(r'^ ^ v^) = 1 — sin^ 29m sin^(1.27Am^L/ifi,) (2) 

where — m^m) is expressed in eV^, L in km and the 

neutrino energy E^, in GeV. Notice the dependence on LjEy in the 
argument of the second sin^ term. 

In the VLI case, the two flavor eigenstates \v^), \vt) and the two 
velocity eigenstates |i/|), [z/g) are connected through the mixing angle 
023 (= 9y) in analogy with mass-induced oscillations: 

Wp) = li^Dcos^la-k |z/^)sin0|3 

Wr) = -\v2) sin6'23 + 1 ^ 3 ) cos0^3 

In this case, the survival probability is 

P{vp = sin^ 29^ sin2(2.54-10^®Aw LE^) (4) 

where Au = ) is the neutrino MAV difference in units of c. 

Notice that neutrino flavor oscillations induced by VLI are character¬ 
ized by an LE,j dependence of the oscillation probability (Eq. 0, to 
be compared with the L/Ei^ behavior of mass-induced oscillations (Eq. 

El- 

When both mass-induced transitions and VLI-induced transitions 
are considered simultaneously, the muon neutrino survival probability 
can be expressed as inmniin] 

P{Vfj, —> Vfj) = 1 — sin^ 20sin^ fl (5) 

where the global mixing angle 0 and the term Q are given by: 

20 = atan(ai/a2) 

_ ( 6 ) 
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with 


ai = 1.27 lAm^ + 2-lO^®Az;sin20„ LE^, e®’*! 

(7) 

02 = 1.27 (Am^ cos29mL/Ei, + 2-lO^®Az;cos20„ LE^) 

Here Am^, L and E^, are expressed, as in Eq. [3 and Eq. 01 in eV^, 
km and GeV, respectively. The additional factor e®®* connects the mass 
and velocity eigenstates, and for the moment it is assumed to be real (77 
= 0 or tt). Note that formulae|21and[nido not depend on the sign of the 
mixing angle and/or on the sign of the Ao and Am^ parameters; this is 
not so in the case of mixed oscillations, where the relative sign between 
the mass-induced and VLI-induced oscillation terms is important. The 
whole domain of variability of the parameters can be accessed with the 
requirements Am^ > 0, 0 < < 7r/2, Au > 0 and — 7r/4 < Oy < tt/A. 

The same formalism also applies to violation of the equivalence 
principle, after substituting Au/2 with the adimensional product [(/jAy; 
Ay is the difference of the coupling constants for neutrinos of different 
types to the gravitational potential m 

As shown in nm, and more recently in im , the most sensitive tests 
of VLI can be made by analysing the high energy tail of atmospheric 
neutrinos at large pathlength values. As an example. Fig. ^ shows the 
energy dependence of the ^ survival probability as a function of 
the neutrino energy, for neutrino mass-induced oscillations alone and 
for both mass and VLI-induced oscillations for Av = 2 x 10“^® and 
different values of sin 29y parameter. Note the large sensitivity for 
large neutrino energies and large mixing angles. Given the very small 
neutrino mass (m^ ^ 1 eV), neutrinos with energies larger than 100 
GeV are extremely relativistic, with Lorentz y factors larger than 10^^. 

MACRO 221 was a multipurpose large area detector (~10000 m^ sr 
acceptance for an isotropic flux) located in the Gran Sasso underground 
Lab, shielded by a minimum rock overburden of 3150 hg/cm^. The 
detector had global dimensions of 76.6 x 12 x 9.3 m^ and used limited 
streamer tubes and scintillation counters to detect muons, z/^’s were 
detected via charged current interactions + N —> ^ + X] upgoing 
muons were identified with the streamer tube system (for tracking) and 
the scintillator system (for time-of-flight measurement). Early results 
concerning atmospheric neutrinos were published in na and in PP for 
the upthroughgoing muon sample and in na for the low energy semi- 
contained and upgoing-stopping muon events. Matter effects in the 
Vsteriie channel were presented in [2] and a global analysis of all 
MACRO neutrino data in |2]. 

In order to analyse the MACRO data in terms of VLI, we used 
a subsample of 300 upthroughgoing muons whose energies were esti¬ 
mated via Multiple Coulomb Scattering in the 7 horizontal rock ab¬ 
sorbers in the lower apparatus in El. The energy estimate was ob¬ 
tained using the streamer tubes in drift mode, which allowed to consid¬ 
erably improve the spatial resolution of the detector (~ 3 mm). The 


3 


overall neutrino energy resolution was of the order of 100%, mainly 
dominated by muon energy losses in the rock below the detector (note 
that {Ef^) ~ 0.4 {E^)). Upgoing muon neutrinos of this sample have 
large zenith angles (> 120°) and the median value of neutrino path- 
lengths is slightly larger than 10000 km. 

Following the analysis in Ref. m, we selected a low and a high 
energy sample by requiring that the reconstructed neutrino energy 
should be < 30 GeV and > 130 GeV. The number of events 
surviving these cuts is Niow = 49 and Nhigh = 58, respectively; their 
median energies, estimated via Monte Garlo, are 13 GeV and 204 GeV 
(assuming mass-induced oscillations). 

The analysis then proceeds by fixing the neutrino mass oscilla¬ 
tion parameters at the values obtained with the global analysis of all 
MACRO neutrino data 12]: Am^ = 0.0023 eV^, sin^26*m = 1- Then, 
we scanned the plane of the two free parameters (Au, 0„) using the 
function 


X 
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A,-aVfg(Au,g.;Am2,g^) \ 
(Ti j 


( 8 ) 


where is the number of events predicted by Monte Carlo, a is 

a constant which normalizes the number of Monte Carlo events to the 
number of observed events and cr^ is the overall error comprehensive of 
statistical and systematic uncertainties. 

We used the Monte Carlo simulation described in ca with different 
neutrino fluxes in input mnniEniEi]. The largest relative difference 
of the extreme values of the MG expected ratio Niow/Nhigh is 13%. 
However, in the evaluation of the systematic error, the main sources 
of uncertainties for this ratio (namely the primary cosmic ray spectral 
index and neutrino cross sections) have been separately estimated and 
their effects added in quadrature (see im for details): in this work, 
we use a conservative 16% theoretical systematic error on the ratio 
Niow/Nhigh- The experimental systematic error on the ratio was esti¬ 
mated to be 6%. In the following, we show the results obtained with 
the computation in m 

The inclusion of the VLI effect does not improve the in any 
point of the (Au, 6^) plane, compared to mass-induced oscillations 
stand-alone, and proper upper limits on VLI parameters were obtained. 
The 90% C.L. limits on Au and 9y, computed with the Feldman and 
Cousins prescription m, are shown by the dashed line in Fig. [3 

The energy cuts described above (the same used in Ref. [T7|L 
were optimized for mass-induced neutrino oscillations. In order to 
maximize the sensitivity of the analysis for VLI induced oscillations, 
we performed a blind analysis, based only on Monte Carlo events, to 
determine the energy cuts which yield the best performances. The 
results of this study suggest the cuts < 28 GeV and £10®° > 142 
GeV; with these cuts the number of events in the real data are 
= 44 events and = 35 events. The limits obtained with this 

selection are shown in Fig. Elby the continuous line. As expected, the 
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limits are now more stringent than for the previous choice. 

In order to understand the dependence of this result with respect 
to the choice of the Am^ parameter, we varied the Am^ values around 
the best-fit point. We found that a variation of Am^ of ±30% moves 
up/down the upper limit of VLI parameters by at most a factor 2. 

Finally, we computed the limit on Av marginalized with respect to 
all the other parameters left free to variate inside the intervals: Am^ 
= Am'^±30%, 9m = 9m ±20%, —7r/4 < 9y < tt/4 and any value of the 
phase 7], We obtained the 90% C.L. upper limit jAul < 3 x 10“^®. 

An independent and complementary analysis was performed on 
a sample of events with a reconstructed neutrino energy 25 GeV < 
< 75 GeV. The number of events satisfying this condition is 106. 
A negative log-likelihood function was built event by event and then 
fitted to the data. We allowed mass-induced oscillation parameters to 
vary inside the MAGRO 90% G.L. region and we left VLI parameters 
free in the whole {Av, 9y) plane. The upper limit on the Av parameter 
resulting from this analysis is slowly varying with Am^ and is of the 
order of « 10“^^. 

In conclusion: we have searched for “exotic” contributions to stan¬ 
dard mass-induced atmospheric neutrino oscillations arising from a 
possible violation of Lorentz invariance. We used a subsample of 
MACRO upthroughgoing muon events for which an energy measure¬ 
ment was made via multiple Coulomb scattering. The inclusion of VLI 
effects does not improve the fit to the data, and we conclude that these 
effects are disfavored even at the sub-dominant level. 

The 90% C.L. limits of VLI parameters are |Az;| < 6 x 10“^^ at 
sin20^ = 0 and |Aw| < 2.5 5 x 10“^® at sin20„ = ±1, see Fig. |21 

In terms of the parameter Av alone (marginalization with respect to 
all the other parameters), the VLI parameter bound is (at 90% C.L.) 
|Au| < 3 X 10-25. 

These results may be reinterpreted in terms of 90% C.L. limits 
of parameters connected with violation of the equivalence principle, 
giving the limit I^Ayj < 1.5 x IQ-^®. 

These limits are comparable or better to those estimated using K2K 
and Super-Kamiokande data |H]. 

We acknowledge the cooperation of the members of the MACRO 
Collaboration. We thank several colleagues for discussions and ad¬ 
vise, in particular B. Barish, P. Bernardini, A. De Rujula, G. L. Fogli, 
S. L. Glashow, P. Lipari and F. Ronga. 
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Av = 2 10 “, sin2e„ = 0 




Av = 210'^®, sin2e„ = ±0.3 




Figure 1: Energy dependence of the survival probability for mass- 

induced oscillations alone (continuous line) and mass-induced -|- VLI oscilla¬ 
tions for Av = 2-10“^^ and sin 29jj = 0, ±0.3, ±0.7 and ±1 (dashed lines for 
positive values, dotted lines for negative values). The neutrino pathlength 
was fixed at L = 10000 km and we assumed Am^ = 0.0023 eV^, 9m = tt/4:. 
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Figure 2: 90% C.L. upper limits on the Lorentz invariance violation param¬ 
eter An versus sin 26y. Standard mass induced oscillations are assumed in 
the two-flavor > Vt approximation, with Am^ = 0.0023 eV^ and 9m = 

7r/4. The dashed line shows the limit obtained with the same selection cri¬ 
teria of Ref. nzi to define the low and high energy samples; the continuous 
line is the final result obtained with the selection criteria optimized for the 
present analysis (see text). 
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